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Abstract. We review results of an analysis of nn interactions in S, P and D waves for two-pion effective 
mass from threshold to about 1.4 GeV. In particular we show a recent improvement of this analysis 
above the KK threshold using more data for phase shifts and including the SO wave inelasticity from 
7T7T — > KK. In addition, we have improved the fit to the /2(1270) resonance and used a more flexible P 
wave parametrization above the KK threshold and included an estimation of the D2 wave inelasticity. 
The better accuracy thus achieved also required a refinement of the Regge analysis above 1.42 GeV. We 
have checked that the tttt scattering amplitudes obtained in this approach satisfy remarkably well forward 
dispersion relations and Roy's equations. 
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1 Introduction 

In a previous analysis [1] a set of fits to different data sets 
on tttt scattering was presented together with a detailed 
description of the mathematical methods used in calcula- 
tions. Forward dispersion relations (FDR) were then used 
in order to test the correctness of the amplitudes thus con- 
structed. Remarkably, it was found that some of the very 
frequently used sets of phase shifts do not satisfy FDR 
below 1 GeV. Thus FDR were shown to give strong con- 
straints to fits which, when used later as a constraint, lead 
to an improved and precise representation of tttt scatter- 
ing amplitudes below, roughly 1 GeV. In the regions from 
about 1 GeV to 1.4 GeV there was still some mismatch 
between the real part of the amplitudes and the results of 
dispersive evaluations in [1] (especially for 7r°7r° scatter- 
ing). 

In a subsequent article [2], in order to improve the 
agreement with the constraints given by FDR, we have re- 
analysed the parametrizations of the 50 above KK thresh- 
old, the DO wave and to a lesser extend the P and D2 
waves. In the SO wave we took into account systematically 
the elasticity data from the tttt — > KK reaction [3-7] , in- 
cluded in the fit more data on phase shifts above the KK 
threshold [3-6] and used more a flexible parametrization 
from 0.932 GeV to 1.4 GeV. In the DO wave we have used 
experimental data from [3,6,8], information on low energy 
parameters (the scattering length and slope) and included 
in the fit the width and mass of the /2(1270) resonance as 
given by the PDG [9]. The result is that, for both SO and 



DO waves, we have obtained more accurate parametriza- 
tions with smaller errors compared to those in the previous 
approach [1] . In the P wave we have exploited a more flex- 
ible parametrization between the KK threshold and 1.42 
GeV and in the D2 wave we have included its estimated 
inelasticity above 1 GeV. 

This more accurate determination of the tttt ampli- 
tudes below 1.42 GeV allowed us to refine the Regge anal- 
ysis that had been used in [1,10]. This has been done 
by removing the degeneracy condition a p (0) — ap>(0) 
which thus modifies slightly the central values of the inter- 
cepts ap(0) (by - 11%) and a P >(0) (by - 4%), but yields 
smaller errors than those in [1] . 

We have found that the tttt amplitudes with the new 
parametrizations of phase shifts and inelasticities in the 
S, P and D waves together with the just discussed small 
changes in a p (Q) and cvp'(0) allow for much better fulfil- 
ment of FDR than in [1] . The biggest improve in % 2 (about 
66%) is achieved for the forward 7r°7r° dispersion relation 
and a smaller one (about 15%) for tt°tt + . In the case of the 
forward dispersion relation for isopin 1 in the i-channel a 
very tiny deterioration has been found (% 2 increased by 
about 26%), which is still acceptable, since, considering 
all FDR together, there is a considerable overall improve- 
ment in their fulfilment. It is worth noting that this has 
been achieved despite the improved data fits have smaller 
errors than in [1]. 

We have also tested Roy equations, which, contrary to 
FDR, incorporate s-t crossing, by calculating the differ- 
ence between the real parts of the input amplitudes and 
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those obtained from Roy's equations. We have found that, 
on average, and up to almost the KK threshold, the devi- 
ation from zero is smaller that 1.05 times the value of the 
errors for the 50 wave, smaller that 1.2 for the S2 wave 
and smaller than 0.65 for the P wave. 



2 S, P and D waves at higher energies but 
below 1.42 GeV 

In this section the main features of the new paramateri- 
zations of S, P and D waves between roughly the KK 
threshold and 1.42 GeV are presented. Details of each 
parametrization can be found in [2] . Since the description 
of the S2 wave was not changed in [2] , any information on 
this wave is available in [1]. 



2.1 The SO wave 

In the present approach we obtain both the tangent of the 
phase shifts tan^fj and the inelasticity T]q above 0.932 GeV 
as functions of K-matrix elements 
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where i,j = 1,2 denote n or K respectively, and we set 
the mass scale /i = 1 GeV. All aj , /3, and 7, are de- 
termined from the fit. Note that M x = 0.9105 ± 0.0070 
GeV s imulates th e left hand cut of the K-matrix located 
at 2y/M\ - ml = 0.952 GeV and the pole at M 2 = 
1.324 ± 0.006 GeV is connected with Sq passing through 
270°. The parametrizations: of [1] (below 0.932 GeV) and 
of [2] (above 0.932 GeV) are matched at 0.932 MeV. In 
the fit all data on phase shifts below and above the KK 
threshold [3-6] have been used simultaneously. For t]q, 
data from irir — > KK have been used together with data 
on 7T7T — * 7T7T [3-7]. The resulting fit yields y 2 /d.o.f = 0.6 
and can be seen in Fig. 1. 



2.2 The DO wave 

For this wave we proceeded by fitting simultaneously be- 
low and above the KK a parametrization 
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with w(s) = ^f^^pZ I, but using different and so pa- 
rameters above and below KK threshold. We also required 
both parametrizations to match at ^fs = 2m k , thus elim- 
inating one parameter. In the present approach, the mass 
of the /2(1270) resonance Mf 2 was fixed to the PDG value 
[9]. The P>i parameters have been obtained for those two 
energy regions from fits to experimental data points [3,6, 
8] together with three other constraints: the width of the 
/2(1270) resonance from [9], plus the scattering length and 
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Fig. 1. Phase shifts and inelasticities of the 50 wave fitted 
using the K-matrix approach of [2] (solid lines). Dotted lines 
for the results of [1]. 



the slope parameter calculated from the Froissart-Gribov 
representation. The resulting \ 2 /d.o.f — 0.65. 

The inelasticity is parametrized in the same way as in 
[1] and fitted to the experimental data of ref.[3,6,8] 
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Results of the fits for phase shifts and inelasticities are 
presented in Fig. 2. 

2.3 The P wave 

In the P wave, above the KK threshold we have used a 
more flexible parametrization than in [1] : 
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where Arj is fixed by the phase shift at 2mx which is ob- 
tained from the fit to the pion form factor [11]. We have 
then fitted data from [6,8] obtaining x 2 /d.o.f '. = 0.6 and 
X 2 /d.o.f . = 1.1 for the phase shifts and inelasticity, re- 
spectively. The results are presented in Fig. 3. 
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Fig. 2. The DO wave phase shifts and inelasticities determined 
in [2] (solid lines) and in [1] (dotted line - only for phase shifts) . 
Dark areas denote the errors, which for the phase shifts have 
just the thickness of the line. 
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Fig. 3. Fits to the P wave phase shifts and inelasticity (solid 
lines). Dark areas show the errors of our results. The dotted 
lines represent results obtained in [1]. 



Fig. 4. Results of the fit to the D2 wave. The solid line denotes 
the fit to the experimental data only and the dashed one the 
fit to the data and FDR [1] . For the data enclosed in the figure 
see references in [2]. 



2.4 The D2 wave 

In the D2 wave we have used one single parametrization 
up to 1.42 GeV with four free parameters 
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where A fixes zero of the phase shift near the 7T7r threshold. 
Since the data on this wave are not accurate we have added 
one more constrain using the scattering length calculated 
from the Froissart-Gribov representation [1]. As a result 
we have obtained the fit presented in Fig. 4. 

The lack of experimental data on inelasticity led us to 
estimate it from a model (see ref. [2]) writing 
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with V§ = 1.05 GeV and e = 0.2 ± 0.2. The inelasticity is 
very small and even negligible below 1.25 GeV. 



3 Regge parametrization 

In the analysis of [1,10] the fits were made with the as- 
sumption of "exact degeneracy" of the intercept param- 
eters a p = otp> for p and fa exchange. In our new ap- 
proach this degeneracy has been lifted. As a consequence, 
there was a very small change in the high energy behaviour 
of scattering amplitudes (especially a little for higher en- 
ergies) but, as can be seen in next section, even such a 
small change could be significant given the level of preci- 
sion achieved in our FDR calculation. The energy depen- 
dence of the new scattering amplitudes after eliminating 
the degeneracy is seen in Fig. 5. 
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Fig. 5. The scattering amplitudes calculated with "exact de- 
generacy" (broken lines) and without this condition (solid 
lines). Dark bands stand for uncertainties. 



4 Implementation of forward dispersion 
relations 

The S, P and D waves presented in Section 2 together 
with the improved Regge description in the previous sec- 
tion have been examined in the same way as in [1], by 
checking the FDR's, but without imposing them as con- 
straints. Thus, in Fig. 6 we present the results from the 
amplitudes in [1] obtained from fits to data. In contrast, 
in Fig. 7 we show results using the improved fits given in 
[2] that we are reviewing here. The Foo, i*b+ and It = 1 
names used in Figs 6 and 7 correspond to the FDR's for 
the 7T°7r , 7r°7r + and t-channel isospin 1 scattering ampli- 
tudes, whose full mathematical expressions can be found 
in [1] and [2]. The word "dispersive" denotes results ob- 
tained from the integrals in the FDR's whereas "direct" 
means the real parts evaluated directly from parametriza- 
tions. 

We provide in Table 1 the FDR's averaged y 2 obtained 
over the range from threshold up to 930 MeV or 1420 MeV. 
Note that the modifications in the S and P waves above 
KK threshold, as well as of the D wave, lead to significant 
improvement of accuracy in the FDR for the 7r 7r scatter- 
ing amplitudes when compared with the previous results 
in [1]. The final decrease of the y 2 f° r tt tt + is also due 
to the influence of the new Regge amplitude. Note that 
in the It = 1 case there is a tiny deterioration despite a 
significant \ 2 decrease due to the Regge part. 



X 



1.4 



-1 - 



-2 - 



— PY direct from data 

— PY dispersive from data 




400 



600 



800 
s 1/2 (MeV) 



1000 



1200 



1400 



Fig. 6. The 7r°7r°, 7r°7r + and t-channel 1 forward dispersion 
relations described in previous analysis [1] . The values of the \ 2 
denote averaged values over all 25 points chosen in the energy 
range from the nir threshold to 1.42 GeV. 



5 Tests using Roy's equations 



We here present an advance of our ongoing analisys where 
we test our new 7T7t scattering amplitudes using Roy's 
equations [12-14]: 
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Table 1. Comparison of averaged x 2 f° r different FDR ob- 
tained in previous analysis [1] and in presented one (new S, r\ 
and new Regge) in two energy ranges. Numbers correspond to 
fits to experimental data only (without constraints from FDR) . 
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Fig. 7. Dispersion relations for new S, P and D waves de- 
scribed in this paper. The \ 2 definition as in Fig. 6 
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with Oq and a 2 , being the SO and S2 scattering lengths, 
KfJ/ (s, s') known kernels and d\(s : s max ) the so called 
driving terms. In our calculations we have chosen s max — 
103m 2 . 

In Fig. 8 we show the real part of the 5*0, S2 and 
P partial waves obtained from Eq.(8), (continuous line, 
called Roy 011 ') versus the real part obtained directly from 
our parametrizations, (dashed line, called called Roy m ). 

The agreement is remarkable, taking into account the 
uncertainties (the dark areas in Fig. 8). Furthermore, the 
agreement is even more impressive, taking into account 
that we have not imposed any constraints from FDR or 
Roy's equations themselves and that the amplitudes come 
just from fits to data (that is why they are labeled "from 
data" in the Figure). Moreover, we use the new S, P and 
D waves described in Section 2 and the Regge model with 
different intercepts a p (0) and otpi (0), and all of them have 
experimental errors even smaller than those of [1]. 



6 Conclusions 

The results reviewed here indicate that the improvement 
in the fits to data in the S and P waves above K K thresh- 
old and the D wave described in Section 2 together with a 
slight improvement in Regge trajectories, allowing for non 
p — f degeneracy, also improves the fulfillment of forward 
dispersion relations. Despite the smaller errors of those 
amplitudes, the averaged x 2 is indeed lower than in pre- 
vious analysis [1]. We have shown that those amplitudes 
fulfill also quite well Roy's equations, and therefore, cross- 
ing symmetry, up to roughly 1 GeV. 

Following, however, the analysis done in [1] one can 
think about a wider implementation of FDR including 
them into the fits together with the already fitted exper- 
imental data. We report briefly on our progress in this 
approach, where, for the moment, we allow for a varia- 
tion of all the amplitude parameters except, the P wave 
above KK threshold and just the a and /3 parameters in 
the K— matrix. Although our results are just preliminary, 
we already noticed significant decreases of the averaged 
X 2 for all three FDR. The preliminary values for Foo de- 
creased from 1.63 to 0.42, for Fq + changed slightly from 
1.44 to 1.48 and for I t = 1 decreased from 1.76 to 0.89. The 
more spectacular improvement, however, has been noticed 
in the Roy equations. Preliminary averaged x 2 values de- 
creased from 1.03 to 0.47 for the 5*0 wave, from 1.18 to 0.53 
for the S2 wave and from 0.65 to 0.02 for the P wave. This 
improvement can be clearly seen when comparing Figs 8 
and 9. Imposing the constrains from Roy equations and 
particularly from FDR, which is much stringent, leads to 
modifications in the S, P and D waves by less than la 
(with the exception for D2 wave where the empirical fit 
changed by ~ 1.3<r) and to negligible modifications in all 
other waves. The resulting uncertainties are also signifi- 
cantly reduced with this approach as can be seen, just for 
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Fig. 8. Differences between real parts of amplitudes calculated 
directly from amplitudes and those from the integral represen- 
tation of Roy's equations. The notation is explained in the 
Section 5. 

Roy equations, in Figs 8 and 9. At present we are finish- 
ing the determination of the final parameters and their 
uncertainties. 



References 

1. J. R. Pelaez, F. J. Yndurain, Phys. Rev. D71, (2005) 
074016. 

2. R. Kamiriski, J. R. Pelaez, F. J. Yndurain, Phys. Rev. D74, 
(2006) 014001 and Erratum Phys. Rev. D74 (2006) 079903. 

3. B. Hyams, et al, Nucl. Phys. B64, (1073) 134. Estabrooks, 
P., and Martin, A. D., Nucl. Physics, B79, (1974) 301. 

4. G. Grayer, et al, Nucl. Phys. B75, (1974) 189. 



Fig. 9. As in Fig. 8 but for data fits constrained by FDR and 
Roy's equations. 

5. R. Kamiriski, L. Lesniak, K. Rybicki, Z. Phys. C74, (1997) 
79 (1997) and Eur. Phys. J. direct C4, (2002) 4. 

6. B. Hyams, et al, Nucl. Phys. B100, (1975) 205. 

7. 7T7T -> KK scattering: W. Wetzel et al, Nucl. Phys. B115, 
(1976) 208; D. Cohen et al, Phys. Rev. D22, (1980) 2595; 
E. Etkin et al, Phys. Rev. D25, (1982) 1786. 

8. S. D. Protopopescu, et al, Phys Rev. D7, (1973) 1279. 

9. Eidelman, S., et al., Phys. Letters B592, 1 (2004). 

10. J. R. Pelaez and F. J. Yndurain, Phys. Rev. D 69, 114001 
(2004). 

11. J. F. de Troconiz, F. J. Yndurain, Phys Rev. D65, (2002) 
093001; Phys Rev. D71, (2005) 073008. 

12. S. M. Roy, Phys. Lett. B 36, 353 (1971). 

13. B. Ananthanarayan, G. Colangelo, J. Gasser and 
H. Leutwyler, Phys. Rept. 353, 207 (2001). 

14. R. Kamiriski, L. Lesniak, B. Loiseau, Phys. Lett. B 551, 
(2003) 241. 



